Performance of the Stereoscopic System of the 
HEGRA Imaging Air Cerenkov Telescopes: 
Monte Carlo Simulations & Observations 



A. Konopelko a , M. Hemberger a , F. Aharonian a , A. Daum a , 
W. Hofmann a , C. K6hler a , H. Krawczynski a , H.J. V61k a , 
A. Akhperjanian b , J. Barrio c ' d , K. Bernl6hr a , H. Bojahr f , 
J. Contreras d , J. Cortina d , T. Deckers 6 , S. DenninghofT c , 
J. Fernandez c ' d , V. Fonseca d , J. Gonzalez d , V. Haustein g , 

G. Heinzelmann g , G. Hermann a , M. HeB a , A. Heusler a , 
H. Hohl f , I. Holl c , D. Horns g , R. Kankanian a ' b , M. Kestel f , 

J. Kettler a , 0. Kirstein e , H. Kornmayer c , D. Kranich c , 

H. Lampeitl a , A. Lindner 8 , E. Lorenz c , N. Magnussen f , 
H. Meyer f , R. Mirzoyan c ' d , H. M611er f , A. Moralejo d , 

L. Padilla d , M. Panter a , D. Petry c ' f , R. Plaga c , 
A. Plyasheshnikov a , J. Prahl 8 , C. Prosch c , G. Piihlhofer a , 
G. Rauterberg 6 , W. Rhode f , A. Rohring 8 , V. Sahakian b , 
M. Samorski 6 , J. Sanchez d , D. Schmele 8 , F. Schroeder f , 
W. Stamm e , B. Wiebel-Sooth f , C. A. Wiedner a , M. Willmer 6 , 
H. Wirth a , (HEGRA collaboration) 

31 Max- Planck- Institut fur Kernphysik, D- 69029 Heidelberg, Germany 

b Yerevan Physics Institute, 375036 Yerevan, Armenia 

c Max- Planck- Institut fur Physik, Fohringer Ring 6, D- 80805 Munich, Germany 

d Facultad de Ciencias Fisicas, Universidad Complutense, E-28040 Madrid, Spain 

e Universitat Kiel, Inst, fur Kernphysik, D-42118 Kiel, Germany 

{ BUGH Wuppertal, Fachbereich Physik, D-^2119 Wuppertal, Germany 

g Universitat Hamburg, II. Inst, fur Experimentalphysik, D-22761 Hamburg, 

Germany 



Preprint submitted to Elsevier Preprint 



1 February 2008 



Abstract 

Based on the Monte Carlo simulations we have studied the performance of the 
HEGRA system of imaging air Cerenkov telescopes (IACTs) in its present config- 
uration of 4 IACTs as well as in its future final configuration of 5 IACTs. Here 
we present the results on the basic characteristics of the IACT system which are 
used in the standard data analysis procedure, i.e., the collection areas, the detection 
rates, the angular resolution, the energy resolution, and the 7/hadron-separation 
efficiency. By comparing several key Monte Carlo predictions with experimental re- 
sults it is possible to check the accuracy of the simulations. The Monte Carlo results 
concerning hadron-nuclear showers are tested with the recorded cosmic ray events 
and the results concerning photon-induced showers are tested with a large data sam- 
ple of 7-rays observed from BL Lac object Mkn 501 during its high flaring activity 
in 1997. Summarizing the simulations and current observations we give the basic 
recommendations of using the instrument and the major values of its sensitivity. 

PACS: 95.55.Ka; 95.55.Vj; 96.40.Pq 

Keywords: Imaging Air Cerenkov Technique; Very High Energy Gamma Ray As- 
tronomy; 



1 Introduction 

The HEGRA collaboration is close to completing an array of five imaging air Cerenkov 
telescopes (IACTs) located on the Roque de los Muchachos, Canary Island La Palma 
(28.8° N, 17.9°). The telescope array, primarily designed for stereoscopic observations of 
the 7-radiation at energies above several hundred GeV [1] is formed by five identical 
IACTs - one at the center, and four others at the corners of a 100 m by 100 m square area. 
The multi-mirror reflector of each system telescope has an area of 8.5 m 2 . Thirty 60 cm 
diameter front aluminized and quartz coated spherical mirrors with focal length of about 
5 m are independently installed on an almost spherical frame of an alt-azimuth mount. 
The FWHM of the point spread function of the reflector is better than 10 arcminutes. 
Each telescope is equipped with a 271-channel camera [2] with a pixel size of about 0.25° 
which results in the telescope's effective field of view of ~ 4.3°. The digitization of the 
PMT pulses is performed with a 120 MHz FADC system. The system trigger demands at 
least two neighboring pixels above the threshold in each of at least two telescopes. The 
present system of four telescopes is taking data in the stereoscopic mode since 1996 [3]. 
The whole telescope array of 5 IACTs is planned to run as a single detector at the end of 
1998. 

The basic concept of the HEGRA IACT array is the stereoscopic approach based on 
the simultaneous detection of air showers by > 2 telescopes, which allows precise re- 
construction of the shower parameters on an event-by-event basis, superior rejection of 
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hadronic showers, and effective suppression of the background light of different origins 
(night sky background, local muons, etc.). The recent observations of the Crab Nebula 
[3] and Mkn 501 [4] by the IACT system strongly support the theoretical expectations 
concerning the features of the stereoscopic imaging technique [5,6]. 

Due to the lack of a calibrated very high energy 7-ray beam, detailed Monte Carlo sim- 
ulations are usually used for the design studies as well as for performance studies of the 
imaging atmospheric Cerenkov experiments. For example, new data analysis methods 
are developed and are tested with Monte Carlo simulations before being applied to real 
data. The measurement of absolute 7-ray flux and energy spectra of the established 7-ray 
sources as well as the determination of upper limits for the quiet objects heavily rely on 
Monte Carlo predictions of the detector performance. In the past, the comparison of the 
characteristics of recorded cosmic ray (CR) events with the characteristics of the Monte 
Carlo simulated hadron-induced air showers used to be the most reliable way to test the 
predictive power of the Monte Carlo simulations (e.g., [7]). Once the telescope response 
to CR-induced air showers is well understood, the Monte Carlo predictions for the 7-ray- 
induced shower can be performed with a high degree of confidence. The situation changed 
dramatically with the observation of the high flaring activity of Mkn 501 in 1997. The 
observations with the HEGRA system of 4 IACTs (~ 110 h) provided a large data base 
of 7-rays (> 30000) with unprecedented signal-to-noise ratio. Several observational key 
characteristics of 7-ray- induced air showers can be measured with small statistical uncer- 
tainties. The agreement of these key characteristics of 7-ray induced air showers in data 
and Monte Carlo simulations substantially strengthened the reliability of the simulations. 

In this paper we discuss the standard data analysis procedure and the results of the 
detailed Monte Carlo simulations for the currently operating system of 4 IACTs as well 
as for the complete HEGRA array of 5 IACTs. Special attention has been paid on the 
proper simulation of the camera and electronics of the Cerenkov telescopes (see Section 
2). Detailed comparisons of the detected cosmic ray and 7-ray-induced air showers with 
simulations have been made to understand the performance of the detector (Section 3). 
The basic characteristics of the HEGRA IACT system were calculated (Section 4). Finally, 
we discuss the resulting sensitivity of the HEGRA IACT array with respect to TeV 7-rays 
(Section 5) and summarize the basic recommendations for the use of the instrument. 



2 Simulations 

In the present calculations the complete Monte Carlo simulation procedure was divided 
into two steps. In the first step an extended library of simulated air showers, induced 
by different primaries, is produced. In the second step simulation of the response of the 
telescope camera pixels is applied to all generated events. The most time consuming step 
is evidently the first one, whereas the second step is relatively fast. This division allows to 
apply the detector simulation procedure several times to the generated showers in order 
to tune the Monte Carlo simulations to the hardware status of the telescopes (e.g. trigger 
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threshold, mirror adjustment etc.) Here we discuss the major features of this two-step 
simulation procedure. 

2. 1 Shower Simulation 

The generation of the air showers has mainly been performed using the ALTAI Monte 
Carlo code [7]. For the electromagnetic cascade this code has implemented an algorithm 
based on the analytical probability distributions of the electron (positron) transport in 
the multiple-scattering segments. This algorithm substantially reduces the computational 
time needed for simulation of a single shower. The proton-nuclei cascade in the atmosphere 
is simulated according to the radial scaling model (RSM) based on accelerator data [8]. 
The air showers induced by the primary nuclei are simulated in a model of independent 
nucleon interactions for the fragmentation at the projectile-nucleus. The fragmentation of 
the colliding nuclei is processed according to the probabilities of different fragmentation 
channels. We studied the influence of the proton-nuclei cascade model on the observable 
shower parameters using the additional simulations with the CORSIKA code [9] which 
has implemented the HDPM model for the simulation of the nucleus-nucleus interactions. 

The shower simulation is carried out at the level of single Cerenkov photons. A certain 
fraction (k ~ 0.2) of the Cerenkov photons of a shower which hit the telescope reflector 
are stored with full information, i.e. the arrival time, the arrival direction, and the impact 
coordinates in the reflector frame. By this means it is possible to apply the complete 
detector simulation procedure to all showers which have been processed in this way. The 
Monte Carlo library contains air showers induced by the primary 7-rays, protons, helium 
and other nuclei belonging to CNO, heavy and very heavy nuclei groups. The primary 
energy of the showers is randomly distributed inside 14 fixed energy bins covering the 
energy range from 100 GeV to 100 TeV. The events are used with weights according 
to some chosen primary spectra. The simulations have been performed for zenith angles 
0°, 30°, 45°. For each type of primary particle and for each zenith angle approximately 
10 5 showers were simulated. The actual setup of HEGRA IACT telescopes [3] was used in 
the simulations. The position of the shower axis in the observation plane was uniformly 
randomized over the area limited by the radius Rq with respect to the central telescope. 
The radius Rq was chosen between 250 and 450 m, increasing with shower energy and 
inclination angle. For the CR air showers the additional randomization over the solid 
angle around the telescope axis with the half opening angle of 3.5° has been introduced 
in order to reproduce the isotropic distribution of the CR events over the camera field of 
view. 

2.2 Detector Simulation 

The detector simulation procedure accounts for all efficiencies involved in the process of 
the Cerenkov light propagation which starts with emission of a Cerenkov photon and 
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ends with the digitization of the PMT signal (see for details [10]). The list of the effects 
which are important in this respect contains: (i) mirror reflectivity, modelled with the ray- 
tracing technique or in phenomenological way using the measured functions of the light 
spot distortion in the camera focal plane; (ii) the light absorption in the plexiglass panel 
covering the camera; (iii) the acceptance of the funnels placed in front of the photomulti- 
pliers (PMTs); (iv) the photon-to-photoelectron conversion inside the PMTs (EMI 9083R) 
taking into account a measured single photoelectron spectrum. The overall efficiency of 
the photon-to-photoelectron conversion is of ~ 0.1. By analogy with the experiment the 
structure of the readout based on the 120 MHz FADC data acquisition and the multiple- 
telescope trigger scheme [11] were implemented in the Monte Carlo simulations. Note, 
that this procedure takes into account the arrival times of the Cerenkov light photons 
hitting the telescope reflector. The basic characteristics of the performance of the tele- 
scope hardware, for instance the number of camera pixels read out for the triggered and 
none triggered telescopes, the single pixel rate and the single pixel trigger rate, the ratio 
between first and second maximum pixel amplitudes in the image etc have been directly 
compared between Monte Carlo and data [10]. 



3 Parameters of the Cosmic Ray Air Showers 

The detection rate of the IACT system is mainly determined by the isotropic flux of 
primary cosmic ray protons and nuclei. The system is triggered if the shower produces 
a sufficiently high number of Cerenkov photons to trigger at least two telescopes. Since 
the number of Cerenkov photons produced by a shower is to first order approximation 
proportional to the energy of the primary particle, the trigger condition determines the 
energy threshold of the IACT system. At larger shower axis distances the Cerenkov photon 
density decreases rapidly. With increasing energy of the primary particles more and more 
showers are able to fulfill the trigger criteria, although they have impact distances far 
away from the telescope system. However, due to the steep primary energy spectrum of 
cosmic rays dJ cr /dE oc E~ 2 - 75 the contribution of the high energy showers (E > 20TeV) 
to the total cosmic ray detection rate is rather small even though they are collected over 
a larger area. 

3.1 Cosmic Ray Detection Rate 

The basic characteristics of a single Cerenkov telescope at hardware level can be calculated 
using the procedure described in [12]. For a system of IACTs the only difference stems 
from the performance of the multi-telescope trigger. For this purpose, first the local trigger 
condition at each telescope as 2nn/271 > qo (the signals in at least two neighboring pixels 
exceed go) has to be fulfiled. The trigger threshold go is measured in photoelectrons. 
The system trigger demands for each individual event coincident trigger signals from 
at least N telescopes (N/5, N = 2, 3, 4, 5). Assuming the energy spectrum and the 
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chemical composition of the primary cosmic rays [13], the cosmic ray detection rate can 
be computed with the following formula: 

n R E 

Rcr = j da2l jR d Rj lm ^P cr (R,E,e ) (!) 



where E is the primary energy of a cosmic ray, R is the impact distance of the shower 
induced by the cosmic ray, 6 is the angle of the shower axis with respect to the telescope 
axis, and P cr (R, E, 9) is the probability of the shower to trigger the telescope system. 

The performance of the multi-telescope trigger of the HEGRA IACT system was discussed 
in detail in [11]. In Table 1 we show the hardware detection rates for the different trigger 
thresholds as measured with the HEGRA telescope system [11] together with the results of 
the Monte Carlo simulations. The measured and computed rates are in a good agreement 
which confirms that the trigger procedure has been modelled correctly. Note that the 
absolute accuracy of the measured rate is about 0.2 Hz. The Monte Carlo predicted rates 
strongly rely on the used chemical composition and fluxes of the cosmic rays. So the 
accuracy of the Monte Carlo predicted rates is better than 20 %. The calculations of CR 
detection rates have been done also for the complete HEGRA system of 5 IACTs (see 
Table 2). One can see that the hardware detection rate for the complete array is expected 
to be approximately ~ 1.4 times larger than for the currently operating 4 IACT system. 

3.2 Cosmic Ray Shower Images 

The standard method to parametrize the Cerenkov light images was primarily introduced 
in [14]. It is based on the second moment analysis of the two-dimensional angular dis- 
tributions of Cerenkov light flashes, sampled with pixels (PMTs) of finite solid angular 
extension [15]. The effective technique (supercuts) to extract the shower image from the 
measured matrix of the pixels amplitude was suggested in [16]. This method provides an 
effective 7/hadron separation and has been extensively used by several single Cerenkov 
telescopes around the world. For the system of IACTs the performance of the method can 
be substantially improved, as will be shown below. 

The distributions of the second-moment image parameters depend crucially on the amount 
of background light per pixel. A slight overestimate or underestimate of the background 
light content dramatically change the distributions. The detailed detector simulation pro- 
cedure is provided to account the exact background content in the camera pixels. In Figure 
1 we show the distributions of the second-moment parameters of the cosmic ray images 
measured in the OFF region of Mrk 501 observations (dark sky region) by the HEGRA 
system of 4 IACT telescopes and the Monte Carlo simulations. It can be seen in Figure 1 
that the Monte Carlo simulated images fit the measured images quite well. 

In general, the distributions of the shape parameters of the Cerenkov light images depend 
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on the model of the development of the proton-nuclei cascade in the atmosphere and also 
on the model of the nucleus-nucleus interaction. To study this effect the series of simula- 
tions have been performed using the ALTAI code [7] and CORSIKA code (with HDPM 
model) [8] in the fragmentation model as well as under assumption of a simple super- 
positon model of nucleus-nucleus interactions. The results show a good agreement of two 
simulation codes, despite of the very different models which are used for the simulations 
of the proton- nuclei component of the air showers. The distributions shown in Figure 1 
have been produced assuming a certain chemical composition of the primary cosmic ray 



4 Imaging of the Gamma-Ray Air Showers 

The hardware detection rate of the cosmic ray air showers dominates by at least two orders 
of magnitude over the detection rate of the 7-rays. To extract the 7-ray signal at sufficient 
confidence level, a special analysis is used to suppress significantly the rate of the cosmic 
ray events. This analysis is based on the application of several software cuts related to the 
orientation and the shape of the Cerenkov light images. Assuming an integral flux and 
an energy spectrum index of a 7-ray source, the detection rates of the 7-ray-induced air 
showers before and after application of the software cuts can be calculated. 



4-1 Collection Areas and Detection Rates 

One of the major advantages of the ground based Cerenkov technique in comparison 
with satellite observations is that the VHE 7-ray-induced air showers can be detected 
at large distances (> 100 m) of the shower core from the telescope. That yields a high 
detection rate of the 7-ray-induced air showers which are distributed over the large area 
of S-f ~ 10 9 cm 2 around the telescope site. The collection area for the 7-ray-induced air 
showers is calculated as 



where P 1 (E,r) is the trigger efficiency of 7-ray-induced air showers of primary energy 
E and impact distance r. The collection area S^(E) for showers of primary energy E, 
is mainly determined by the effective area of the telescope reflector S p h. e . = S m ■ Xph.e., 
where S m is the total mirror area and Xph.e. is the efficiency of the photon-to-photoelectron 
conversion of the camera channels. Given a fixed mirror area, a maximum collection area 
is achieved by reducing the trigger threshold of the telescopes. The latter is limited at the 
lower end by the fluctuations of the background light in each camera pixel. 



[13]. 



00 




(2) 
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In Figure 2 the collection areas for the complete HEGRA system of IACTs are shown 
for the conventional trigger criteria. The strong increase of collection area in the energy 
range < 1 TeV changes to a logarithmic growth at higher energy. The behaviour of the 
collection areas is determined by the shape of the lateral distribution of the Cerenkov 
light pool at the observation level. The density of Cerenkov light photons at plateau 
for impact distances up to 125 m (for air showers observed at angles close to zenith) is 
roughly proportional to the primary shower energy whereas beyond 125 m the density 
of the Cerenkov light density decreases rapidly. In observations at large zenith angles 
(> 30 degree) the collection area decreases at low energies (E < 3 TeV) but increases at 
higher energies (E>5 TeV). The reason for that is that air showers at the larger zenith 
angles develop higher in the atmosphere. This is clue to the fact that the same amount of 
the Cerenkov light (neglecting the increase in absorption) is produced by a shower, but 
is scattered over the larger area at the observation level which provides decreasing the 
density of Cerenkov photons. Thus low energy air showers at large zenith angles cannot 
trigger the telescopes but at the same time the high energy air showers have much larger 
collection area due to the large size of the Cerenkov light pool at observation level. 

The HEGRA system of 5 imaging air Cerenkov telescopes has been designed for effective 
observation of 7-rays with the primary energy of several hundred GeV in stereoscopic 
mode with telescopes of relatively small mirror area, 8.5 m 2 . The system trigger based on 
the simultaneous detection of the shower images in several telescopes (at least 2) forces 
down the trigger threshold and consequently the energy threshold of the IACT system. 
Usually, the energy threshold is determined as the energy at which the detection rate of 
observed 7-ray showers reaches its maximum. For convenience we use in the following rate 
calculations a 7-ray spectrum according to: 

dJ^/dE = A ■ E~ a \ J 7 (> 1 TeV) = KT 11 cnrV 1 . (3) 

For a certain spectrum index a 7 the 7-ray detection rate is calculated as 

00 00 



where (dR 1 /dE) [Hz TeV -1 ] is the differential 7-ray detection rate. 

Under the assumption of a differential index, a 7 , of the 7-ray energy spectrum one can 
calculate the differential 7-ray detection rate (dR^/dE) (see Figure 3). The peak of the 
differential detection rate slightly shifts to higher energies with increase of the trigger 
multiplicity because low energy events cannot effectively trigger several telescopes. The 
integral detection rates of the 7-ray-induced air showers for different trigger multiplicities 
are presented in Table 3. Note that in the case of a steep 7-ray spectrum (e.g. a 7 = 3.0) 
for the system trigger 2/5 the 7-ray detection rate is much higher as compared with the 
operation for the 3/5 trigger mode. It is seen from Figure 3 that for the trigger conditions 
2/5, 2nn/271> q ph.e. the differential detection rate peaks at an energy of ~ 500 GeV 
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which is identified as the energy threshold of the instrument. The energy threshold for 
observations at large zenith angles significantly increases (see Figure 4). At the same 
time, the rate of high energy events detected at large zenith angles could even exceed 
the corresponding rate in observations at the nominal zenith angles due to the significant 
increase of the collection area with increasing zenith angle. 

Our Monte Carlo studies show that 7-ray air showers detected at large impact distance 
from the telescopes cause some difficulties for a reliable selection of 7-ray events. For 
plane-parallel 7-ray flux the Cerenkov light image shifts to the camera edge for large 
impact distances. These images are partially cut by the camera edge and cannot be used 
for an accurate reconstruction of the shower parameters (orientation of the shower axis 
in space, shower core location etc.). Thus for better evaluation of the energy spectrum 
it is useful to set a restriction on the reconstructed impact radius from the center of the 
system. This restriction influences mainly the collection areas at high energies. Above a 
certain energy the effective collection area is then determined simply by the geometrical 
area around the IACT system. The restriction on the impact distance for r < 200 m in 
the case of steep spectrum does not significantly change the detection rate. For the case of 
a flat energy spectrum, especially at observations at large zenith angles, it is an advantage 
to collect 7-ray-induced air showers at large distances to the center of the IACT system. 



4-2 Reconstruction of Shower Arrival Direction 

The simultaneous observations of the air showers with > 2 imaging air Cerenkov telescopes 
offers the possibility to reconstruct the orientation of the shower axis with respect to the 
telescope axis [18]. The general approach is based on the superposition of the several 
images in one common focal plane in order to derive the intersection point of the major 
axis of the ellipsoid-like images. This intersection point determines the shower direction 
[6]. If the Cerenkov telescopes are directed towards the object, the reconstructed source 
position in the camera field of view for the 7-ray-induced air showers has to be in the center 
of the camera focal plane. The currently operating HEGRA system of IACTs performs 
so-called wobble mode observations. The position of the source in the camera focal plane is 
offset by 0.5° from the camera center (on declination) and consequently rotates depending 
on the azimuth angle. This approach gives the possibility to perform continuous ON source 
observations, whereas the OFF region can be chosen in 1° offset from the source position 
[3]. In present simulations the Cerenkov light images were shifted by 0.5° from the center 
of the focal plane with the correlated randomization over the azimuth. 

The difference between the true and reconstructed position of the 7-ray source in the 
camera field of view, O, is a measure of the angular resolution of the system of IACTs. 
The distributions of O 2 , both for the simulated and observed 7-ray showers from Mrk 501 
are shown in Figure 5. One can see from Figure 5 both distributions match and both show 
a prominent peak around the source position. Our Monte Carlo studies show that the tail 
of the distribution at large O 2 is due to the air showers with the core positions close to the 
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line connecting two or three telescopes. In such a case, the images in the telescopes are 
almost parallel to each other and the reconstruction procedure leads to significant error 
in the evaluated shower direction because of a small intersection angle. Note that in the 
reconstruction procedure we require at least 3 telescopes have to be triggered (in addition 
we require also a minimum number of photoelectrons in the image of 40 ph.e.). 

The angular resolution of the system of IACTs can be characterized quantitatively by the 
acceptance of the 7-ray-induced air showers, K^ ir , after the application of the fixed angular 
cut on © 2 . In Table 4 the data on the acceptance, K^ ir , of the 7-rays for three angular cuts 
@ 2 < 0.03, 0.05 and 0.1 [deg 2 ] are shown for the simulations at different zenith angles 0, 
30, 45 degree. In general, for observations at large zenith angles the shower is far from 
the observer and the Cerenkov light images detected at large zenith angles have a smaller 
angular size, they are closer to the camera center and show almost circular shape. These 
changes in the image topology lead consequently to larger errors in the reconstruction of 
the shower direction. The comparison of the angular resolution for two different trigger 
multiplicities, 2/5 and 3/5, show some increase in the 7-ray acceptance applying the same 
angular cuts for the higher trigger multiplicity - 3/5. However, because of the energy 
dependence of the angular resolution two different system triggers have to be considered 
as complementary in observations of 7-ray sources with very different spectral features. 

The angular resolution noticeably depends on the impact distance of the shower core from 
the center to the I ACT system. For impact distances r < 125 m the angular resolution 
slightly improves with increasing the energy of the 7-ray showers, because the images on 
average contain more light and the image orientation is better determined. The angular 
resolution, 50 (one standard deviation of the Gaussian distribution on 9), for ~1 TeV 
7-ray showers is of 0.11° and 0.09° at ~ 20 and 200 m, respectively. Beyond ~ 120 m the 
angular resolution decreases at higher energies. This can be explained by the distortion 
of the images by the edge of the limited camera field of view. The impact distances for 
high energy 7-rays correspond to large shifts of the images from the center of the camera 
focal plane (the high energy air showers occur deeper in the atmosphere). For ~ 10 TeV 
7-ray showers detected at impact distances of 200 m the angular resolution is 0.14°. 

The advanced angular resolution of the HEGRA system of IACTs is a very effective tool 
for suppression of the isotropic cosmic ray background. The acceptances of cosmic ray air 
showers after application of an angular cut, @ 2 < Oq, evaluated from the data taken with 
the currently operating HEGRA system of IACTs (OFF sample of Mrk 501) as well as 
from the Monte Carlo simulations are shown in Table 5. It is seen from Table 5 that after 
applying the angular cut of ©q = 0.03 [deg 2 ] the cosmic ray background rejection is as 
high as ~ 200. One can see also from Table 5 that Monte Carlo simulations reproduce 
quite well the measured contamination of the cosmic ray air showers after application of 
an angular cut. 
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4-3 Localization of Shower Core 



The position of the shower core at the observation level can be measured by the system 
of IACTs for a single individual event. Then the impact distances from the shower core 
to the system telescopes can be evaluated. The reconstruction algorithm is based on the 
orientation of the Cerenkov light images in several telescopes which have been triggered. 
We use pure geometric reconstruction which does not relate to the image shape. The 
accuracy of the shower core reconstruction is limited by the errors in the determination 
of the image orientation. As discussed before, the change of the Cerenkov light image 
topology leads to an increase of the error in the core position with increasing zenith 
angle. The accuracies of shower core reconstruction for the different primary energy and 
impact distance are summarized in Table 6. Observing a 7-ray source at zenith angles 
below 45 degrees and restricting the impact distances from the telescope to the shower 
core within 200 m the average accuracy is < 20 m. 

The reconstructed impact distance is used for calculating the image shape parameters 
scaled on an impact distance and image amplitude. These parameters are applied for the 
cosmic ray rejection. For that, the accuracy of 20 m is quite sufficient because the shape 
of the Cerenkov light images does not significantly change within 20 m. Furthermore, the 
value of the reconstructed impact distance is needed for evaluation of the shower energy. 
The calculations show that even with an accuracy of the shower core localization of around 
20 m the energy resolution for 7-ray-induced air showers is better than 20%. 



4-4 Measurement of Shower Energy 

The procedure of energy reconstruction for the 7-ray-induced air showers observed with 
a single imaging air Cerenkov telescope appears to be quite complicated due to the lack 
of a direct measurements of the distance from the telescope to the shower core. Several 
indirect methods have been invented in order to estimate the impact distance using the 
centroid shift in the camera field of view as well as the image shape [19,20]. 

For the system of IACTs the measurement of the impact distance is straightforward and is 
not related to the image shape. That improves the accuracy in the energy reconstruction 
compared with a single telescope. The algorithms of the energy reconstruction for each 
individual event as well as for the spectrum evaluation with the system of IACTs have been 
discussed in [4,5,21]. This method of energy reconstruction was investigated by means of 
Monte Carlo simulations for the HEGRA system of IACTs. 

If one can measure the distance from the shower core to the telescope (ri,i = 1,N, 
where N is the number of triggered telescopes) the primary energy of air showers can be 
evaluated using the inverse function of the image size (amplitude) on shower energy and 
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impact distance 



(5) 



where Si is the image size (total number of photoelectrons in the image) and 9 is a zenith 
angle. The final energy estimate can be constructed by incorporating several images in 
a number of telescopes, E^i = 1,N, as E = J2i"WiEi, where Wi is an energy an core 
distance dependent weight (J2 w i — !)• Note that the energy resolution depends on the 
accuracy of determing the core distance and is strongly influenced by fluctuations in the 
image amplitude. 

In Table 7 we show the estimates of the energy resolution in the different primary energy 
ranges as well as for the simulations at different zenith angles. Note that in a case of 
wobble mode observations the images of high energy 7-ray-induced air showers E > 10 TeV 
are very often cut by the camera edge. This leads to a distortion in the impact distance 
reconstruction and consequently in the reconstruction of the shower energy. To remove this 
effect the large impact distances r > 200 m are usually excluded from the consideration. 
For inclined showers (at large zenith angles) the Cerenkov light images are closer to the 
center of the field of view because these showers develop in atmosphere very far from the 
telescope and the corresponding images shrink to the camera center. Thus the problem 
of the camera edge appears to be less significant for inclined showers in the energy range 
at least up to 20 TeV. Data shown in Table 7 demonstrate that over the whole energy 
range for the 7-ray-induced air showers in observations at zenith angles up to 45 degrees 
the estimated energy resolution for the HEGRA system of IACTs is around 20%. 

4-5 Rejection of CR background using the image shape 

Due to the difference in the nature of 7-ray- and proton/nuclei-induced air showers the 
corresponding Cerenkov light images are also very different in shape [22] . In the standard 
second-moment approach these differences can be determined by using such image shape 
parameters as Width, Length etc. The 7-ray-induced air showers have on average, images 
of smaller angular size. The selection of the 7-rays can be done by applying the standard 
parameters cuts: Width < wo, Length < l etc where wo, Iq are the boundaries which limit 
the domain of the most of the 7-ray-induced showers (> 50%). 

When the minimum amplitude for the pixel used for the image parametrization is fixed, 
the parameters of the image shape would increase with the total number of photoelectrons 
in the image. For the high energy air showers more pixels are involved in the image 
parametrization procedure that is why the angular size of the image increases with the 
shower energy. Application of fixed image shape cuts saves the most of the low energy 
7-rays (which are close to the energy threshold) and reduces significantly the content of 
high energy 7-rays. For the 7-ray spectrum evaluation that is undesirable because it will 
decrease the statistics of the observed high energy 7-rays. To avoid this problem the energy 
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dependent cuts have to be used. In addition, for a fixed primary energy of 7-ray-induced 
air showers the angular size of image also depends on the distance from the telescope to 
the shower core due to the decrease in image amplitude with impact distance. Thus for 
a single telescope using the restriction on the position of the image in a certain range 
from the camera center, one can keep only events at core distances r < 120 m. The radial 
dependence of the angular size of an image is very small in this particular range of impact 
distances. However, such a restriction significantly reduces the number of detected 7-rays. 
For the system of IACTs both the energy and radial dependence of the image size can 
be accounted in detail, because the IACTs system measures the shower core position and 
consequently the energy for each individual event. 

The images from the 7-ray-induced air showers can be sorted into several bins on the 
measured distance from the shower core, Ar«, i — 1, n, and image size, Alog(Sj),j = 1, m. 
Then the averaged image shape parameters are calculated for these particular bins, < 
w >ij, < I >ij. For an individual event the shower core position and impact distance from 
the shower core to the triggered telescopes can be reconstructed. Instead of the usual 
Width (w) of the image the scaled Width [17] (w) is calculated for each telescope and 
after that the mean scaled Width parameter is defined for this event: 

N 

<w>= l/N Y^w k / <w >% (6) 
k=i 

where N is the number of triggered telescopes. The rejection of cosmic ray background 
events is performed by applying a cut on the mean scaled Width, < w ><< w > . 
Measured distributions of the mean scaled Width for 7-ray-induced air showers observed 
at different zenith angles < 10, 30, 45 degrees are shown in Figure 6. The distributions 
have a peak at 1.0 and are very narrow because after the scaling procedure the RMS of 
the distributions is determined only by the pure image fluctuations and does not depend 
any more on the image amplitude and shower impact distance. The multifold cut on the 
image Width is replaced by single parameter mean scaled Width and works effectively 
against the cosmic ray background events. For a single telescope (i.e. one view of the 
shower) the probability, that cosmic ray shower gives an image beeing very similar to the 
7-ray-induced shower, could be as high as 0.1. In observations by several telescopes (in 
different projections) this probability is reduced down to 10~ 2 (see Table 8). Note that 
in addition to mean scaled Width cut another shape parameters, such as Length, could 
be used. However, the simulations show that already after the mean scaled Width cut 
the cosmic ray background rejection is already very strong and an application of another 
shape cut does not improve the signal-to-noise ratio but simply reduces the amount of the 
7-rays. The acceptances of the cosmic ray showers after the application of a mean scaled 
Width cut can be compared for Monte Carlo simulations and observed cosmic ray events. 
The data shown in Table 8 demonstrates a good agreement between the simulations and 
data. 

There are two possible strategies for the choice of cuts. The first one is to apply the so- 
called strong cuts (< w >= 1.0) which dramatically suppress the CR background but at 
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the same time also noticeably reduce the content of 7-rays (see Table 8). This approach 
works best when searching for the 7-ray signal from source candidates on a short time 
scale. For spectrum studies it is more effective to apply the loose cuts (< w >= 1.4) which 
save most of the observed 7-rays at high energies and do not show strong energy dependent 
efficiencies (opposite to the tight cuts). It is important to note that the efficiency of 
the cosmic ray rejection improves approximately by a factor of 2 for 3-telescope events 
compared to 2-fold coincidence events while the decrease of the 7-ray acceptance is only 
of ~ 20%. Thus use of 3-fold events seems to be preferable from the point of view of 
the cosmic ray rejection using the image shape and orientation as well as for an accurate 
impact distance and energy reconstruction. However, in the specific case of the 7-ray 
fluxes characterized by very steep energy spectra the analysis of 2-fold coincidences could 
be applied in order to increase the statistics of 7-rays. 



5 Sensitivity to 7-ray fluxes 

Using the calculated characteristics of the IACT system performance, one can estimate 
the sensitivity of the instrument in 7-ray observations. For two complementary approaches 
of the cosmic ray background rejection based on the application of tight and loose cuts the 
acceptances of the 7-rays and cosmic rays are shown in Table 9. For the integral flux of 
7-rays taken at the level of J 7 (> 1 TeV) = 1CT 11 cm~ 2 s _1 the expected number of detected 
7-rays per one hour is about 28 and 56 for the tight and loose cuts, respectively. For the 
tight cuts the cosmic rays are highly suppressed and the expected cosmic ray rate per hour 
is around 3 particles. Thus, the stereoscopic observations of 7-ray sources of an integral flux 
J 7 (> 1 TeV) > 1CT 11 cm _2 s _1 with the HEGRA array of IACTs are essentially background 
free (rate of the detected 7-rays exceeds the rate of the background event by more than a 
factor of 10). The expected sensitivity of the 5 IACT system for the 7-ray flux stated above 
after application of the loose cuts is about 5 "sigma-per-one-hour" . The 5a detection of 
7-ray fluxes from point sources at the level of J 7 (> 1 TeV) = 10 -12 cm~ 2 s _1 is possible in 
20 hours of observations. Note, that after application of the tight cuts the amount of the 
detected cosmic electrons (~ 0.4 particle per hour) becomes comparable to the cosmic 
ray contamination and further suppression of the background is limited by the electron 
content. 



6 Summary 

The HEGRA system of imaging air Cerenkov telescopes is the first instrument operating 
in the stereoscopic observation mode. The HEGRA system of 5 telescope with relatively 
small mirror area of 8.5 m 2 gives a rather low energy threshold of ~ 500 GeV. The 
large camera field of view (~ 4.3 degree) allows to perform the observations of the point 
7-ray sources in the wobble mode, taking at the same time the ON and OFF events, 
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and increase the available observation time by a factor of two. The use of several images 
gives the angular resolution better than 0.1 degree and yields the cosmic ray rejection, 
using the shower orientation, up to 200 times. The geometrical reconstruction of the 
shower impact position with a good accuracy (< 20 m) improves the energy resolution 
and makes possible to apply the image shape cuts which are independent on a shower 
energy. The energy reconstruction procedure for the telescope system is straightforward 
and is not related to the image shape. The energy resolution of the HEGRA system in 
the dynamic energy range of 0.5 GeV - 30 TeV is better than 20 %. In order to avoid the 
uncertainties in the evaluation of the energy spectrum of the detected 7-rays the effective 
collection area can be taken simply as the geometrical area around the center of the 
telescope system using the restriction on the reconstructed impact parameter (e.g., 200 
m). Simulateneos registration of a several Cerenkov light images from the individual air 
shower makes possible to apply the correlated analysis of the image shape (mean scaled 
Width) and substantially improve the cosmic ray rejection up to 100. Note that the data 
analysis based on three concidence view appears to be preferable for the better angular, 
energy resolution as well as gives better cosmic ray rejection using the image shape. The 
three concidence events are optimum for the stereo imaging of the TeV 7-ray air showers. 
In the search mode the tight software cuts on shape and orientation show the better 
performance approaching the case of a background free detection of the 7-ray sources. 
However for the energy spectrum studies one can use the loose cuts providing the high 
7-ray statistics. 

The HEGRA system of IACTs could be considered as a successful prototype for the 
future low energy (100 GeV) arrays such as HESS and VERITAS. In general the current 
experience of the HEGRA operation can be used for the design of the forthcoming arrays 
of IACTs. 
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Fig. 1. Distributions of image parameters of the recorded cosmic ray air showers and of Monte 
Carlo simulated cosmic ray air showers at the zenith angles up to 30 degrees. The hatched 
histograms correspond to the data, the points are for the Monte Carlo simulations. 
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Fig. 2. Collection areas of the 7-ray- induced air showers simulated at zenith for the complete 
array of 5 HEGRA IACTs. Data correspond to different trigger multiplicities N = 2, 3, 4 and 5 
(N is a number of triggered telescope for an individual event). Local trigger for each telescope 
is taken as 2nn/271> 10 ph.e. 
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Fig. 3. The differential detection rates of 7-ray-induced air showers simulated at 0° zenith angle 
and detected by the system of 5 HEGRA IACTs operating with different local trigger condition: 
2/5, 2nn/271 > qo, <7o= 8 (1); 10 (2); 14 (3) ph. e. The assumed 7-ray spectrum is a power-law 
with the spectral index of 2.5 and with a flux normalization chosen according to Eqn. 4. 
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Fig. 4. The collection areas (a) and differential detection rates (b) for the 7-ray-induced air 
showers simulated at zenith angles 0°, 30°, 45° and system trigger: 3/5, 2/271>10 ph.e. The 
assumed 7- ray spectrum is a power-law with the spectral index of 2.5 and with a flux normal- 
ization chosen according to Eqn. 4. 
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Fig. 5. Distribution of reconstructed arrival directions of the 7-ray-induced air showers. is an 
angular distance of the reconstructed position of the source in the camera focal plane from the 
true source position. The histogram corresponds to the Monte Carlo simulations. The distribu- 
tion for 7-rays detected from Mrk 501 is indicated by black dots. 
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Fig. 6. Distributions of Width (left column) and mean scaled Width (right column) parame- 
ter of the 7-ray-induced air showers detected by the HEGRA IACT array (dots) at different 
zenith angles ~0, 30, 45 degrees. Curves and d^stograms show the results of the Monte Carlo 
simulations. 



Table 1 

The cosmic ray detection rates, R cr [Hz], for the currently operating HEGRA system of 4 IACTs 
in observations at the zenith angles up to 30 degree. The data correspond to the different multi- 
telescope triggers: N/4 (N = 2, 3, 4) in the case of trigger 2nn/271 > qo for each individual 
telescope. The measured rates are taken from [11] (data are interpolated for the particular trigger 
thresholds). 
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Table 2 

The cosmic ray detection rates, R cr [Hz], for the complete system of 5 HEGRA IACTs calculated 
at the zenith angles up to 30 degree for different multi-telescope triggers, N/5, and for two 
trigger thresholds go = 8, 10 ph.e. 
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Table 3 

The detection rates of the 7-ray- induced air showers, i? 7 [photons/hour], for the complete 
HEGRA system of 5 IACTs for the different multi-telescope triggers. The local trigger for each 
individual telescope is set as 2nn/271 > qo, qo= 10 ph.e. The results correspond to a differential 
7-ray spectrum in form of a power law with different spectrum index 07. 
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Table 4 

The Monte Carlo acceptance, , of the 7-ray-induced air showers simulated at different zenith 
angles, 6, after application of the angular cut Q 2 < Gq [deg 2 ]. The calculations have been done for 
the system trigger of 2/5 and 3/5 and the local telescope trigger 2nn/271>10 ph.e. Additionally, 
the restriction on the distance of the shower core from the central telescope < 200 m was applied. 
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Table 5 

Acceptances of cosmic ray air showers, K^ r , after application of the angular cut, G 2 < Gq- Data 
and Monte Carlo (MC) simulations correspond to zenith angles up to 30 degree, system trigger 
3/5 and the local telescope trigger 2nn/271>10 ph.e. 
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Table 6 

Accuracy in localization of shower core of the 7-ray-induced air showers, 5r (one standard 
deviation of the Gaussian approximation), simulated at different zenith angles, 6, and impact 
distances from the center of IACT array, r. The impact distance ranges are numbered as 0-50 
m (1); 50-100 m (2); 100-150 m (3); 150-200 m (4); 200-250 m (5). 
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Table 7 

Energy resolution, AE/E, for 7-ray-induced air showers simulated at the zenith angles of 0, 30, 
45 degrees. The maximum impact distance is 200 m. Data correspond to a system trigger of 3/5, 
and a local telescope trigger 2nn/271>10 ph.e. 
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Table 8 

The acceptance of the 7-rays, re 7 , the cosmic ray background contamination, K cr , after application 
of a cut on mean scaled Width, < wq >. The Monte Carlo simulations and data correspond to 
a zenith angle range up to 30 degree. 
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Table 9 

The expected rates of 7-r ay-induced and cosmic ray induced air showers, simulated at zenith 
angles up to 30 degree, before and after application of the angular and shape cuts. The estimated 
signal-to-noise ratio for the 7-ray source with the intensity of J 7 (> ITeV) = 1CT 11 cm -2 ^ 1 
(power-law spectrum index was taken as a 7 = 2.5). Data corresponds to the system trigger 3/5, 
the local telescope trigger 2nn/271>10 ph.e. and 200 m maximum distance from the shower 
axis. 

Tight cuts (W <1.0, @l < 0.05 [deg 2 ]J: 

i Ri, Hz nf r n s j hape Ri, hour- 1 
7 1.7 -10- 2 0.86 0.53 28 
CR 7.8 0.009 0.01 3 

Loose cuts (W <1.3, B 2 , < 0.1 [deg 2 ]/: 

i Ri, Hz 4 ir K f ape Ri, hour' 1 
7 1.7 -10 -2 0.92 0.99 56 
CR 7.8 0.018 0.15 76 
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